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SUMMARY 

Transient heat conduction in a slab of finite thickness is analyzed so as to include the 
effect of position- dependent heat generation. One face of the slab is a convecting surface, 
and the other is considered to be insulated. Temperature distributions are obtained for 
the cases of constant, linear, exponential, and cosinusoidal heat- production rates. The 
solutions are obtained by the integral transform method and are in terms of Fourier se- 
ries. Numerical values of these series were calculated and are presented graphically, 
in dimensionless form, so that specific cases can be easily computed. The method of 
solution is outlined in detail so that temperature distributions can be obtained for other 
position- dependent heat- generation rates. 


INTRODUCTION 

Heat generation in solids occurs in the fields of electrical, chemical, and nuclear 
engineering. Closed-form analytical solutions for the resulting conduction temperature 
distributions are generally limited to simple- shaped bodies with constant heat- production 
rates. 

Transient heat conduction in a slab with one insulated surface and constant heat pro- 
duction is treated in references 1 to 3. Solutions, with numerical results, are given in 
references 1 and 3 for the case of a slab initially at zero temperature with the uninsu- 
lated surface maintained at zero temperature. The solution to the same problem except 
for a slab with a constant initial temperature is given in reference 2. Also in refer- 
ence 3, charts are presented for the surface temperature variations when the uninsulated 
surface is convecting and the slab is initially at the temperature of the surrounding me- 
dium. 

The conduction problem in the slab initially at zero temperature, with one insulated 


and one zero-temperature surface, is developed formally in reference 1 for an arbitrary 
position- dependent heat- generation rate. The resulting temperature distribution is in the 
form of an integral, and specific generation functions are not treated. 

In the field of nuclear engineering, the heat- generation function is usually not a con- 
stant but is a function of position relative to the reactor core (refs. 4 and 5). The de- 
tailed form of the generation- rate function depends on factors such as the core configura- 
tion. The variation is often cosinusoidal or exponential, and in some cases it can be ap- 
proximated by a linear variation. 

Consideration is given herein to transient heat conduction with position-dependent 
heat generation in a finite slab that has one insulated and one convecting surface. The 
slab is initially at a constant temperature. The problem is developed for heat generation 
with an unspecified variation in the direction of the slab thickness. The cases of constant, 
linear, exponential, and cosinusoidal heat- production variations are considered in detail. 


ANALYSIS 

Formulation of Problem 

Consider the transient conduction of heat in a homogeneous finite slab with constant 
properties k, p, and c (fig. 1). (Symbols are defined in appendix A. ) Initially the slab 
is at the constant temperature T^. For t > 0, heat is generated in the slab at a rate 
G(X), which depends on position in the direction normal to the bounding planes. Convec- 
tion heat transfer occurs at the face X = 0 between the body and the surrounding medium, 
which is at the constant bulk temperature T^. The convective heat-transfer coefficient 
is considered constant, and the plane X = L is thermally insulated. Then the governing 
equation, the initial condition, and the surface conditions are the following (ref. 1): 

= 0 <x<l <d 

at 9X 2 pc 

T = T. t = 0 for all X (2) 

k H = h(T - T h ) at X = 0, t > 0 (3) 

3X d 


— = 0 at X = L, t > 0 (4) 

ax 
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Solutions 


The problem defined by equations (1) to (4) can be solved by the method of integral 
transform. A modified Fourier cosine transform is applied to the space variable. This 
procedure yields an ordinary differential equation in time and an initial condition for the 
transformed temperature. Let 


x = L - X 

e = t - T b 


Then equations (1) to (4) can be rewritten as follows: 


90 - a d2 ° + 0 < x < L 

(la) 

9t 3x 2 P c 


= T i " ^ = 0 i t = 0 for all x 

(2a) 

— = 0 at x = L, t > 0 

(3a) 

3x k 


^=0 at x = 0, t > 0 

(4a) 


3x 

Applying the modified Fourier cosine transform (see appendix B) to equations (la) and 
(2a) and noting that Cj = Cg = 0 give 



dt l 2 





pc 


(5) 


rss _ 

0 = 0. t = 0 


( 6 ) 


where 
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0 . cos dx = sin X„ 

1 L X n 
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Ax 

G = / G(x) cos dx 
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( 8 ) 


and A n are the roots of 


A_ tan A_ = Bi 
n n 


(B3) 


The solution to equation (5) with the initial condition (eq. (6)) is 


-X*?Fo / -X?Fo\ p. 2 
" + |l-e n ' GL 


kX 
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( 9 ) 


o 

where Fo = at/L is the Fourier number. 
Inverting equation (9) gives 


e = 



(A n + sin A n cos A n ) 


0 i e 


-xJfo 
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- X n Fo \ GL 


kX 


n 


cos — A 
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( 10 ) 
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which in terms of the original variables is 


T - T b 

T i” T b 



(A n + sin A n cos A n ) 
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e n sin A n + 1 - e n 9J± 
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n 


x cos 


[(-W 


(10a) 


Thus, for a specific position- dependent heat- generation rate, the transform of the gen- 
eration function can be evaluated using equation (8). The transient temperature distribu- 
tion is then given by equation (10a). The cases of constant, linear, exponential, and co- 
sinusoidal heat- generation rates are treated herein. 

Constant heat- generation rate. - For a constant heat- generation rate (fig. 2), 


4 



t > 0 


( 11 ) 


G(X) = G q 


Then, from equation (8), 



( 12 ) 


Substituting equation (12) in equation (10a) gives the following temperature distribution 
for constant heat production: 


— 4 * i + - 


Gq l 


T, - T u 1 k(T t - T b ) 




(13) 


The position- and time- dependent functions and herein referred to as the tran- 
sient temperature and uniform temperature functions, respectively, in equation (13) are 



2 

where Fo = ort/L is the Fourier number and Bi = hL/k is the Biot number. The 
values of X n in equations (14) and (15) are the roots of equation (B3), namely, 

X n tan X n = Bi (B3) 


Linearly varying heat- generation rate . - When the heat- generation rate varies lin- 
early (fig. 3), the following equation results: 

G(X) = G Q + mX t > 0 (16) 
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where m is the constant (positive or negative) slope of the generation rate. The trans- 
formed generation function (eq. (8)) for this case is 


G = 


f [< G o * 

t/0 


Ax G n L T 2 

mL) - mxjcos -5- dx = sin X n + (1 - cos Aj 


n 


n 


n' 


(17) 


and the temperature distribution (eq. (10a)) becomes 


T - T h G n L T 3 

b = ^ + 0 xj. + — 

T. - T b k(Tj - T b ) 2 k(T. - T b ) ‘ 


(18) 


The linear temperature function \p^ (eq. (18)) is defined as 


-A?Fo 


1 - e n 1(1 - cos A n ) 
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(A„ + sin X cos A J 
v n n n' 
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(19) 
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Exponentially varying heat- generation rate . - The heat- generation function for ex- 
ponential variation (fig. 4) is 


G(X) = Gpe'^^/^ t>0 


( 20 ) 


where ju. is a. constant attenuation coefficient (positive or negative). Substituting equa- 
tion (20) in equation (8) gives 


G = 


/ 


G 0 e 


-/i[l-(x/L)] x n x . 
v / /J cos -ii- dx = 


G 0 L 

2 .2 
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Me‘ M ) (21) 


Then, from equation (10a), the temperature distribution is 
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( 22 ) 
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where 



(23) 


is the exponential temperature function, and A n are the roots of equation (B3). 

Cosinusoidally varying heat- generation rate . - The cosinusoidal heat- generation 
rates indicated in figure 5 are 


G(X) = G m cos — X 
M 2D 


(24) 


The positive or negative 6 in equation (24) is the displacement of the convecting surface 
from the plane of maximum heat production. 

With cosinusoidal heat generation, the transformed generation function is given by 


G = 


/ 




cos — f(6 + L) - xlcos 
2D 
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dx 


(25) 


Evaluating the integral gives 


G = 


G M L 
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sin(d + A) - d sin A cos A n - cos A sin Aj 


when X^ * d 2 (25a) 
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where 


d-li 

2 D 


A 71 6 

A = 


2 D 


Then the temperature distribution with cosinusoidal heat generation is 


v, g m l 

b = +— 


T - T 

T i -T b "’' 1 T k<T 1 -T b ) 


The cosinusoidal temperature function in equation (26) is 
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- d sin A cos X n - X n cos A sin Xjcos 
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and 


when X 2 * d 2 (27) 


* 5 = dcos(d + A) + sindcosA L _ e -d 2 Fo\ CQS 1 _ xV 1 
d 2 (d + sin d cos d) V / l\ W J 


when X 2 = d 2 (28) 


RESULTS AND DISCUSSION 

Numerical values of the dimensionless time- and position- dependent ifs functions 
were calculated according to equations (14), (15), (19), (23), and (27). The first six 
roots of the transcendental equation (eq. (B3)) were used for all calculations. These 
roots were taken from reference 1. The results are presented graphically in figures 6 
to 10. 

Individual charts are given for values of dimensionless position X/L from 0 to 1.0 
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in increments of 0. 2. The range of Fourier numbers covered is 0 to 1. 4. The Biot num- 
bers included are 0, 0. 2, 0. 5, 1. 0, 2. 0, 4. 0, 10. 0, and °°. A Biot number of 0 corre- 
sponds to an insulated surface. Similarly, the infinite Biot number results apply for the 
case of a prescribed constant surface temperature T^. 

The transient temperature function xp^ is given in figure 6. This function is the di- 
mensionless transient temperature distribution with no heat production that is commonly 
presented in the literature (e. g. , refs. 1 and 6 to 9). The uniform temperature function 
i ^2 is given in figure 7. The charts given in reference 3, namely, the temperature- 
distribution chart for the slab with prescribed surface temperature and the surface- 
temperature chart for the slab with a convecting surface, are special cases of figure 7. 
The numerical results are in agreement throughout the range of Fourier number (0 to 0. 2) 
presented in reference 3. The linear temperature function xp^ is presented in figure 8. 
These results apply for positive or negative heat- production- rate slopes. The exponential 
temperature function xp^ depends on the attenuation coefficient p as well as on X/L, 

Fo, and Bi. Values of xp ^ are given in figure 9 for p values of 0. 5, 1. 0, 2. 0, 4. 0, 
-0.5, -1.0, -2.0, and -4. 0. Similar to the preceding case, the cosinusoidal temperature 
function xp g depends on two additional parameters, d and A. Values of xp^ for the rel- 
ative slab thickness d = 7 t/ 6 and the displacement of the uninsulated face A = -n/3, 

-7t/ 6, 0, and 77/6 are given in figure 10. 

For small slab thicknesses, the Fourier number is large even for moderate values of 
time. Also, for small L, the Biot number will be low except for large values of the con- 
vective coefficient. For cases of small thickness, however, the variation of heat gener- 
ation with position will be small, and uniform production can be assumed. Values of the 
transient temperature function xp ^ and the uniform temperature function xp^ for large 
Fourier numbers (0 to 5. 0) and small Biot numbers (0 to 4. 0) are given in figures 11 
and 12, respectively. 

The range of the values presented is indicated in the following table for a few specific 
materials and slab thicknesses: 


Material 

Slab thickness 

Temperature 

function 

Maximum 

time, 

min 

Minimum convective 
heat-transfer coefficient, 
h, for steady state 

in. 

cm 

Btu/ (hr ) (ft 2 ) (°R) 

J/(hr)(m 2 )(°K) 

Stainless steel 

4 

10.2 

All 

62 

96 

1. 96x1 0 6 


2 

5. 1 

Transient 

55 

48 

.98 


2 

5. 1 

Uniform 

55 

48 

. 98 

Aluminum 

12 

30. 5 

All 

25 

466 

9. 4 xlO 6 


4 

10. 2 

All 

3 

1400 

28.6 


2 

5. 1 

Transient 

2. 5 

700 

14. 3 


2 

5. 1 

Uniform 

2. 5 

700 

14. 3 
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Figures 6 to 12 should be adequate for most engineering calculations. For values of the 
parameters not included in these figures, the temperature functions can be computed di- 
rectly from the series forms presented herein. 

As an example of the effect of heat- generation position dependence on the conduction 
temperature distributions, consider the case of an aluminum rod 1 foot (0. 305 m) long 
under the following conditions: 

h = 240 Btu/(hr) (ft 2 ) (°R) ; 4. 9xl0 6 J/(hr)(m 2 )(°K) 

T b = 200° R; 111 0 K 

T. = 200° R; 111°K 

Gq = 200 000 Btu/(hr)(ft 3 ); 74.6X10 8 J/(hr)(m 3 ) 

k = 120 Btu/ (hr) (ft) (°R) ; 0.746X10 6 J/(hr)(m)(°K) 

a = 3.33 ft 2 /hr; 0. 310 m 2 /hr 
o 

Then, B- = hL/k = 2.0 and Fo = at/L = 0. 0555 t, where t is in minutes. 

The case of exponential heat generation with ju. = 1. 0 is compared with the case of 
uniform heat generation with a constant value of Gq. From equations (22) and (13), 

G n L 2 

T = T b +— 5 — j// 4 = 200 + 1670 (°R); 111 + 927 i^ 4 (°K) for the exponential case 

k 

G L 2 

T = T fe + -° — xf/ = 200 + 1670 (°R); 111 + 927 \j/ 2 (°K) for the uniform case 

k 

The values of ^ 4 and are obtained from figures 9 and 7, respectively. 

The time-temperature history at the insulated end of the rod is shown in figure 13. 
The calculated temperatures are considerably higher under the condition of uniform gen- 
eration. Also, the time predicted to reach a specified temperature (say 860° R or 478° K) 
is much less for the uniform- generation case. 
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The temperature distributions at t = 15 minutes are shown in figure 14. At all posi- 
tions in the rod, the predicted temperature values with uniform generation are higher 
than values calculated for exponential heat production. Figure 14 also shows that temper- 
ature gradients under the condition of uniform generation are larger than with exponential 
heat generation. 


CONCLUDING REMARKS 

The transient conduction of heat in a solid bounded by two parallel planes was ana- 
lyzed so as to include the effect of position-dependent heat production. It was assumed 
that the solid is homoegeneous with constant physical properties and constant initial tem- 
perature. The slab with one convecting surface and one insulated surface was considered. 
Heat- production rates varying constantly, linearly, exponentially, and cosinusoidally in 
the direction normal to the bounding planes were treated in detail. 

The solutions in terms of Fourier series were obtained by the integral transform 
method. The temperature distributions depend on position, time (Fourier number), con- 
vection coefficient (Biot number), and parameters that reflect the specific form of the 
heat- generation rate. Numerical values of the temperature-distribution series were cal- 
culated and are presented graphically in dimensionless form. The range of variables 
was so selected that specific cases can be easily computed. 

The solution procedure is presented for a heat- generation rate with an unspecified 
variation in the direction of the slab thickness. Thus, by following the outlined proce- 
dure, transient temperature distributions can be determined for specific position- 
dependent heat- generation rates other than those treated in detail in this report. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 20, 1967. 
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APPENDIX A 


SYMBOLS 


Bi 

Biot number, hL/k 

A 

dimensionless displacement, 

C-p Cg 

constants 


Or/2) (6/D) 

c 

specific heat 

5 

displacement of convecting face 

D 

quarter wavelength for cosine 
generation 

e 

for cosinusoidal heat generation 
temperature difference, T - 

d 

dimensionless slab thickness, 


roots of eq. (B3) 


Or/2) (L/D) 

M 

attentuation coefficient for expo- 

e 

exponential base 


nential heat- generation rate 

Fo 

2 

Fourier number, art/L 

P 

density 

f 

modified Fourier cosine trans- 

*1 

transient temperature function 

G 

form defined by eq. (Bl) 
volumetric heat- generation rate 

*2 

defined by eq. (14) 
uniform temperature function de- 

h 

convective heat-transfer coeffi- 
cient 

^3 

fined by eq. (15) 
linear temperature function de- 

k 

L 

thermal conductivity 
slab thickness 

^4 

fined by eq. (19) 

exponential temperature function 
defined by eq. (23) 

m 

slope of linearly varying heat- 

'I'*. 

cosinusoidal temperature function 

T 

generation rate 
temperature 

u 

defined by eqs. (27) and (28) 

Subscripts: 

T b 

bulk temperature of surroundings 

i 

initial 

t 

time 

M 

n 

maximum 
1, 2, 3, . . . 

X 

coordinate in direction of slab 

X 

a 

thickness 
coordinate L-X 
thermal diffusivity, k/pc 

0 

value at X = 0 
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APPENDIX B 


MODIFIED FOURIER COSINE TRANSFORM 

The modified Fourier cosine transform used to solve the system of equations (la) 
to (4a) is defined (refs. 10 and 11) as 



(Bl) 


where the function f(x) must satisfy the Dirichlet conditions (i. e. , bounded with a finite 
number of maxima and minima and a finite number of points of discontinuity). The bound- 
ary conditions of the problem are 


— = C. at x = 0 
dx 1 

► 

= -- f + C 9 at x = L 
dx k > 


In equation (Bl) A n are the roots of 


A„ tan A„ 
n n 



The transform of the Laplacian of f turns out to be 






Cj + C 2 cos A n 


The inversion formula is 



(B2) 


(B3) 


(B4) 


(B5) 
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Convecting Insulated 

surface surface 



Figure 1. - Finite slab. 



Coordinate in direction of slab thickness, X 

Figure 2. - Slab with constant heat- 
generation rate. 



Slope of linearly varying 
heat -gene ration rate, m 

Positive 

Negative 


Coordinate in direction of slab thickness, X 
Figure 3. - Slab with linearly varying heat -gene rat ion rate. 



Attenuation coef- 
ficient for ex- 
ponential heat- 
generation 
rate, p 

Positive 

Negative 


Coordinate in direction of slab thickness, X 
Figure 4. - Slab with exponentially varying heat -gene ration rate. 
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Volumetric heat-generation rate, G(X> 



Coordinate in direction of slab thickness, X 

(a) Positive displacement of convecting face. <b) Negative displacement of convecting face. 


Figure 5. - Cosin usoidally varying heat-generation rates where D 
generation, G M is maximum volumetric heat-generation rate, L 
placement of convecting face for cosinusoidal heat generation. 


is quarter wavelength for cosine 
is slab thickness, and 6 is dis- 
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Transient temperature function, ^ 




(b) Dimensionless position, 0.2. 



(d) Dimensionless position, 0.6. 



(e) Dimensionless position, 0.8. (f) Dimensionless position, 1.0. 

Figure 6. - Response of transient temperature function for slab subjected to sudden change in environ- 
mental temperature. 
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Fourier number, crt/L 2 

(e) Dimensionless position, 0,8. (f> Dimensionless position, 1.0 

Figure 7. - Response of uniform temperature function ^ f° r slab with uniform internal heat generation. 









(a) Dimensionless position, 0 


(b) Dimensionless position, 0.2. 



• 2 .4 .6 .8 1.0 1.2 1.4 0 .2 .4 

Fourier number, at/L 2 

(e) Dimensionless position, 0.8. 


.8 1.0 1.2 1.4 


<f) Dimensionless position, 1.0. 

Figure 8. - Response of linear temperature function ^ for slab with heat generation varying linearly in direction of slab 
thickness. 
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Exponential temperature function, ip 4 





(a) Attenuation coefficient, 0.5. 

Figure 9. - Response of exponential temperature function ip 4 for slab with heat generation varying exponentially in direction 
of slab thickness. 
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Exponential temperature function, 






Exponential temperature function, 



Fourier number, ort/L 2 


(e) Attenuation coefficient, -0.5. 
Figure 9. -Continued. 










0 .2 .4 .6 .8 1.2 1.4 1.6 0 .2 .4 .6 .8 1.2 1.4 1.6 

Fourier number, at/L 2 

(a) Dimensionless slab thickness, 7r/6; dimensionless displacement, -7r/3. 


Figure 10. - Response of cosinusoidal temperature function for slab with heat generation varying cosinusoidally in di- 
rection of slab thickness. 
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Cosinusoidal temperature function, ifj 5 




1.2 1.4 0 .2 

Fourier number, crt/L 2 



(c) Dimensionless slab thickness, 7T/6, dimensionless displacement, 0. 
Figure 10. - Continued. 


Cosinusoidal temperature function, $ 5 




(d) Dimensionless slab thickness, tt/6; dimensionless displacement, tt/6. 
Figure 10. - Concluded. 


32 



asssS 


Hi!8i§™alBiSilSSSliM 

laK fflia HWHKSBH HMHWagM 

lgg|IKpS 33 gi|ggSSPBB 3 H 3 | 

iBHHiiiltesgaBBjBaggagiiisggI 


inl'M 


(a) Dimensionless position, 


maun 

liawi— 

IlijgHB gg aii^Wiliiiigilil; 

ISHPSSRBSHd 


IISKBSSNffl^BESS 


israsp 


SKf'IMil 


gaggsgagwaaii^asl 


(b) Dimensionless position, 0.2. 


iHilklKl 

min 


SSsgj 


iBBIia iaHHfiSBaHHH BBiaa— 
iHB&lllBBBB^SBBBBBBMBBhW 


(c) Dimensionless position, 0.4. 


NHSI 


a BiUB BBgB BBBBHBBBBBBM 

■RutbbbbbbbbbbbbbsbbSI 


IBMBB B BBB IBaiB BgB BBHgiBiy 
lillBB B BBlM BBBBaBB3B§i 


(d) Dimensionless position, 0.6. 
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Fourier number, at/L c 


(e) Dimensionless position, 0.8. 


(f) Dimensionless position, 1.0. 


Figure 11. - Response of transient temperature function tp-^ for slab subjected to sudden change in environmental tern 
perature; high Fourier number and low Biot number. 







Fourier number, at IV 


(e) Dimensionless position, 0.8. (f) Dimensionless position, 1.0. 

Figure 12. - Response of uniform temperature function ip ^ f° r slab with uniform internal heat generation; high 
Fourier number and low Biot number. 














